Nonlinear parametric instability in double-well lattices 



Jaroslav Riedel 1 , Boris A. Malomed 2 , and Eva Majerm'kova 1 
1 Department of Theoretical Physics, Palacky University, 
Tf. 17. Listopadu 50, 77207 Olomouc, Czech Republic 
2 Department of Interdisciplinary Studies, Faculty of Engineering, 
Tel Aviv University, Tel Aviv 69978, Israel 

A possibility of a nonlinear resonant instability of uniform oscillations in dynamical lattices with 
harmonic intersite coupling and onsite nonlinearity is predicted. Numerical simulations of a lattice 
with a double- well onsite anharmonic potential confirm the existence of the nonlinear instability with 
an anomalous value of the corresponding power index, ~ 1.57, which is intermediate between the 
values 1 and 2 characterizing the linear and nonlinear (quadratic) instabilities. The anomalous power 
index may be a result of a competition between the resonant quadratic instability and nonresonant 
linear instabilities. The observed instability triggers transition of the lattice into a chaotic dynamical 
state. 



Dynamical lattices with onsite nonlinearity and har- 
monic intersite coupling constitute a vast class of models 
which have numerous physical application and are an obp 
ject of great interest in their own right, see, e.g., Ref.EJ. 
Among these models, the ones with a double-well on- 
site potential [given, e.g., by the expression ( |l2| ) below] 
have special importance, as they directly apply to the 
description of structural transitions in dielectrics, semi- 
conductors, superconductors, and optical lattices (see 
recent worka3 and references therein), and find other 
applicationsB. 

The simplest dynamical state in lattices represents spa- 
tially homogeneous oscillations. This state in conserva- 
tive lattice models is sometimes stable, and sometimes it 
is subject to linear modulational instabilities initiating a 
transition to nontrivial dynamicsQ. An objective of the 
present work is to demonstrate analytically, and verify by 
direct simulations, that homogeneous oscillatory states 
in lattices with the double-well anharmonicity may be 
subject to a specific nonlinear instability, which triggers 
transition of the lattice into a chaotic dynamical state. 

The nonlinear parametric instability, which is a sub- 
ject of this work, is inherently related to the phonon an- 
harmonism in the lattice: the instability is caused by a 
resonance involving the uniform oscillations of the lattice 
and a second harmonic of the phonon mode. A possibility 
of a nonlinear instability of an intrinsic localized morfdll 
in the lattice due to the phonon anharmonism was first 
considered in Ref.B. An opposite, and more common, 
type of the resonance, namely, between a strictly linear 
phonon mode and a higher harmonic of an intrinsic lo- 
calized mode is well known to give rise to a slow decay of 
the localized mode into phonons (see recent worksa and 
references therein). 

The general form of the lattice equation of motion is 

«n + / (««) u n = u n +i + u n -i - 2u„, (1) 

where u n are real dynamical variables on the lattice, the 
overdot stands for d/dt, f (u n ) is a polynomial function 
accounting for the onsite nonlinearity (in fact, nonpoly- 
nomial functions can be considered too), and the right- 



hand side of the equation accounts for the intersite har- 
monic coupling. The linearized version of Eq. (Q) gives 
rise to phonon modes 

u n = Asm (kn — ujt) (2) 

with an arbitrary infinitesimal amplitude A and the dis- 
persion relation uj — 2 |sin(fc/2)| inside the phonon band, 
w < 2 (by definition, the frequencies are positive). 

A homogeneous oscillatory state Uo(t) is a time- 
periodic solution to the equation 

+ / (U ) U = (3) 

with a fundamental frequency Q. The linear stability 
of the homogeneous state is determined by a linearized 
equation for small perturbations 6u ni which is produced 
by the substitution of u n = Uo(t) + 8u n into Eq. (|IJ): 

Su n + [/ ' (U (t)) U (t) + f (U (t))} Su n 
= Su n+ i + 5u n -i - 2Su n . (4) 

In the mean-field approximation, one may describe 
phonon modes of the type (^) on top of the homogeneous 
oscillations, replacing the coefficient in front of Su n on 
the left-hand side of Eq. (0) by its time-average value 
wg = (f(U (t))U {t) + f(U (t))). The corresponding 
dispersion relation for the phonon modes acquires a gap 
ujq, so that 

uj 2 = u% + 4sin 2 (/c/2), (5) 
which gives rise to the phonon band 

< ^ 2 < 4 + w o ( 6 ) 

(if loq < 0, the homogeneous oscillations are immediately 
unstable). 

Beyond this simple approximation, the Fourier de- 
composition of the coefficient / ' (Uo(t)) Uo(t) + f (Uo(t)) 
in Eq. (Q) gives rise to parametrically driven terms 
~ cos(milt) ■ Su n with all integer values of m. The linear 
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parametric drive resonates with a perturbation frequency 
w, i.e., it may give rise to a resonant linear instability, un- 
der the condition mil — u — to, or 

u = u$*> = (m/2)0. (7) 

If any resonant frequency (to/ 2) Q gets into the phonon 
band (||), the homogeneous oscillatory state is expected 
to be modulationally unstable, otherwise the resonant 
linear instability does not take place. 

In the latter case, it makes sense to seek for nonlinear 
parametric instabilities, the simplest of which may be 
generated by a cubic term, or any higher-order one, in 
the onsite nonlinearity. Indeed, the cubic term generates 
a nonlinear correction ~ U~o(t) ■ (5u n ) 2 to Eq. M), which 
may be regarded as a parametric drive that can give rise 
to a nonlinear parametric resonance under the condition 
mQ — 2lu = u>, or 

u = c^° nlin) = (m/3)n, (8) 

where to is an arbitrary integer different from a multiple 
of 3, cf. Eq. (0) (if to is a multiple of 3, the linear 
parametric resonance takes place at the same frequency, 
so that the nonlinear resonance is insignificant). 

Of course, this description has a very approximate na- 
ture for two reasons. First, the phonon band ([|) was 
defined in the framework of the mean-field approxima- 
tion, hence one cannot be sure in the accuracy of the 
predictions based on the comparison of the resonant fre- 
quencies with this band. Second, the full lattice model 
(j]) may give rise to other instabilities, which are not re- 
lated to the parametric resonance. Therefore, the above 
consideration should only be considered as a qualitative 
clue, and an actual possibility of dynamical regimes dom- 
inated by the nonlinear resonance must be checked by 
direct simulations. 

Continuing the consideration, we note that, if none of 
the linear-resonance frequencies (0) gets into the renor- 
malized band but a nonlinear-resonance frequency 
(||) can be found inside the band, an evolution equa- 
tion for the amplitude of the corresponding resonant- 
perturbation mode, 8u n = A(t) cos {u> res t)-v n , with some 
spatial profile v n (it may be, for instance, the above- 
mentioned localized intrinsic mode), has a general form 

dA/dt = CA 2 , (9) 

where C is a constant which depends on a particular 
form of Eq. (Q) and the homogeneous solution Uq (t) ; cf. 
a similar equation governing the nonlinear instability of 
the so-called embedded solitonsn. A solution to Eq. (f|) is 

A = A / (1 - CA t) , (10) 

where A Q is the initial value of the perturbation ampli- 
tude. A drastic difference of the perturbation growth law 
( |l0| ) from the exponential growth in the case of the lin- 
ear instability is that the nonlinear instability is initially 



growing much slower than an exponential, and a char- 
acteristic time scale of the growth, ~ 1/ (CA ), depends 
on the initial perturbation Aq, while in the case of the 
exponential growth it is a fixed constant. However, the 
nonlinear instability is self-accelerating, and, as a mani- 
festation of that, Eq. (|Io| ) formally predicts a singularity 
at t = 1/ (CAq). In reality, of course, the singularity may 
not occur, as the above approximation, taking into regard 
the first nonlinear correction to Eq. (Q), becomes irrele- 
vant if A(t) is too large. A natural conjecture, that will 
be corroborated by direct simulations below, is that the 
nonlinear instability leads to a chaotic dynamical state. 

It is relevant to mention that, although nonlinear in- 
stabilities are less common than the usual linear insta- 
bility, they occur and play an important role in many 
physical problems, as diverse as optical solitons in_media 
with competing quadratic and cubic nonlinearitiesElj-Bose 
gases, plasma turbulence, contact lines in flows, etc.E2l. In 
this work, we will check the possibility of the nonlinear 
instability of the homogeneous oscillations in the lattice 
model (||) with 

/(«„) = -u 2 n + vui, v > 0, (11) 

which corresponds to the double-well onsite anharmonic 
potential, 

V(u n ) = -u 4 n /A + vu 6 n /6. (12) 

In this case, Eq. (||) can be solved in terms of elliptic 
functions, but an explicit result is very cumbersome. 

As a typical example, we take homogeneous oscilla- 
tions produced by Eq. (^) with v = 0.01 and initial 

conditions Uq(0) = 1 and [To (0) = 0. The variable Uq 
then performs strongly anharmonic oscillations between 
the values (U ) min = 1 and (£/o) max = 12.247 at the 
fundamental frequency = 1.694, and the gap in the 
renormalized phonon spectrum (^) is calculated to be 
ujq = 3. 518, so that the renormalized band (|^) is, in the 
mean-field approximation, 

3.518 < lu < 4.047. (13) 

Then, it is straightforward to check that all the linear- 
resonance frequencies (Q) do not get into this band (the 
band as whole fits between the linear resonant frequencies 
3.392 and 4.240, which correspond to to = 4 and to = 5). 
On the other hand, the nonlinear-resonance frequency (|^) 

corresponding to to = 7 is Wre S onlm ^ = 3. 9573, which lies 
inside the band ([l3]) [all the other frequencies given by 
Eq. ( |l3| ) are located outside the band]. 

To directly test the instability, small perturbations of 
the form 

6u n (0) = A Q cos(2TT Po n/N), (14) 

where N is the net number of sites in the lattice and po 
is an integer, were added to the homogeneous oscillatory 
state. The lattice equations of motion were solved for 
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N = 1000 and periodic boundary conditions by means 
of the eighth-order explicit Runge-Kutta scheme with a 
stepsize control such that the time step was dynamically 
changed within the range 0.05 - 0.3. It was checked that 
the relative (per site) error at each step did not exceed 

io- 10 . 

The simulations were performed for the perturbations 
( |l4| ) with po taking values in the interval 1 < po < 30. In 
all the cases considered, results were quite similar. Here, 
we demonstrate a typical example with po = 20. Long- 
time evolution initiated by the small perturbation Jl4| ) 
with Ao = 0.05 is displayed is Fig. 1 in the form of a 
set of plots showing the temporal development of several 
components in the Fourier transform of the lattice field, 
which are defined as follows: 

JV 

U p {t) = [2 /N) Unit) exp (2inpn/N) , p ^ 0; 

n=l 
N 

U (t) = (l/N)Y J U n (t). (15) 



It is obvious that the small perturbation triggers a tran- 
sition of the lattice into a chaotic state. Fully developed 
chaos, i.e., a state in which all the lattice modes are in- 
volved into the chaotic motion, is attained at t w 22, 
when the phonon mode with p = po + 1 gets chaoti- 
cally excited too, see Fig. 1. To further illustrate the 
transition to chaos, in Fig. 2 we additionally show in 
detail, on the logarithmic scale, the growth of the ampli- 
tude \U Po+ i(t)\. More detailed studies of the established 
chaotic state may be of interest in their own right, but 
this problem is beyond the scope of the present work. 

As concerns the nonlinear character of the instability, 
a crucial issue is the growth of the perturbation at the 
initial stage. It is necessary to check whether it is indeed 
essentially different from the familiar exponential law, 
being, instead, close to the Eq. (|l0|). To this end, in 
Fig. 3 we display the best fit of the time evolution of 
the numerically computed Fourier amplitude \U Po (t)\ to 
a function 



A fit (i) = A • (1 - yty 



(16) 



where the parameters are found to be A = 0.041, 7 
0.560 and a = 1.750. 



Comparison of these results with Eq. (10) shows a dif- 
ference in the (most essential) power parameter a. Note 
that the expression ( fl6| ) with the empirically found value 
a = 1.750 formally corresponds to a solution to the non- 
linear evolution equation dA/dt = CA@, with an anoma- 
lous value of the power index, (3 = 1 + or 1 » 1.5714, 
that should be compared to Eq. L ^J), valid in case 
of the ordinary nonlinear instabilitytJ. This anomalous 
value is sort of intermediate between (3 — 1 and (3 = 2, 
which are expected for the linear and nonlinear insta- 
bilities, respectively. This result may suggest that, in 
fact, in the present model we have a competition between 



the resonant nonlinear instability, qualitatively consid- 
ered above, and linear instabilities against nonresonant 
perturbations, which were not taken into regard in the 
above consideration. While an accurate analysis of the 
full linear stability problem of the homogeneous oscilla- 
tions is a technically complex problem, that we do not 
aim to consider here, Fig. 3 clearly shows that the reso- 
nant nonlinear instability dominates in the growth of the 
perturbations. 

In conclusion, we have proposed a possibility of a non- 
linear resonant instability of homogeneous oscillations 
in harmonically coupled nonlinear lattices, which is ex- 
pected to play a dominant role, provided that no resonant 
frequency accounting for the linear parametric resonant 
instability gets into the renormalized phonon band, while 
a frequency that gives rise to a quadratic parametric res- 
onance is found in the band. Numerical simulations of 
the lattice with a double-well onsite anharmonic poten- 
tial confirm the existence of nonlinear instability with an 
anomalous value of the power index ss 1.57, which is in- 
termediate between the values 1 and 2, characteristic of 
the linear and nonlinear instabilities. The onset of the 
nonlinear instability triggers transition of the lattice into 
a chaotic dynamical state. 
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FIGURE CAPTIONS 



Fig. 1. The time dependence for selected Fourier am- 



shown for p = 0, p = po, p = 2po, and p = pq + 1, where 
Pa = 20. The lattice size is N = 1000 (with periodic 
boundary conditions), and the initial amplitude of the 
perturbation is Aq = 0.05. 

Fig. 2. Details of the evolution of the Fourier ampli- 
tude |£/ po +i(t)|, shown on the logarithmic scale. 

Fig. 3. Fitting the time dependence of the amplitude 
\U Po (t)\ to the function @ with A = 0.0411, 7 = 0.560 
and a = 1.750. Diamonds stand for numerical data, 
and stars (which almost completely overlap with the di- 
amonds) show the closest values provided by the fitting 
function. 
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